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W)| Abstract 

We investigate f(R, T) gravity models (R is the curvature scalar and T is the trace of the stress- 
energy tensor of ordinary matter) that are able to reproduce the four known types of future finite-time 

"xl" ' 

singularities. We choose a suitable expression for the Hubble parameter in order to realise the cosmic 

acceleration and we introduce two parameters, a and H B , which characterise each type of singularity. 

-^ ' We address conformal anomaly and we observe that it cannot remove the sudden singularity or the 

type IV one, but, for some values of a, the big rip and the type III singularity may be avoided. We 

also find that, even without taking into account conformal anomaly, the big rip and the type III 



singularity may be removed thanks to the presence of the T contribution of the f(R, T) theory. 
Pacs numbers: 04.50.Kd, 95.35.+d, 95.36. +x, 98.80.Qc 

1 Introduction 

Observation indicates that the current expansion of the universe is accelerating [TJ [5]. There are 
two large branches of cosmology attempting to explain this phenomenon [3]-|12j: i) Introducing dark 
energy in the framework of general relativity is the first; ii) The other is investigating modified versions 
of the gravitational theory. In this paper we follow the latter approach. The most widely renown 
modified theory is /(ij)-gravity, in which the action is determined by an arbitrary function f(R) of 
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the Ricci scalar R [5]-[8]. Another interesting modified theory of gravity is /(G)-gravity, where G = 
R 2 — AR^R^ + R^ pa R^ vp,J is the Gauss-Bonnet invariant (R^ u and R^ vpiJ are the Ricci tensor and the 
Riemann tensor respectively). 

In the framework of general relativity, it is well known that a phantom phase usually ends up in the 
type I (big rip) finite-time future singularity (FTFS) [T3J [14]. It has also been demonstrated, e.g. in 
[TS]-[T7], that the dark energy equation of state parameter may show all four possible types of FTFS 
and that there is no qualitative difference between dark energy and modified gravity. Thus, /(i?)-gravity 
dark energy models also may bring the universe evolution to all four possible FTFS [EJ-J22]. Another 
interesting aspect of f(R) modified gravity is that it may provide models that remove FTFS by adding, 
for example, a i? 2 -term [TH1 EM [HI [23] ■ The /(G)-gravity also is a class of modified theories that may 
lead to finite-time future singularities [HI [24] ■ 

In the present paper, we approach /(i?, T)-gravity, where T is the trace of the energy-momentum 
tensor of ordinary matter, and investigate models that may lead to the four FTFS. This type of modified 
theory of gravity has been first developed by Harko et al. [25] , who derive the gravitational field equations 
in the metric formalism, as well as the equations of motion of test particles, which follow from the 
covariant divergence of the stress-energy tensor. They also analysed the Newtonian limit of the equations 
of motion, and provide a constraint on the magnitude of the extra-acceleration by investigating the 
perihelion precession of Mercury. In [26 , f(R, T) models have been constructed describing the transition 
from the matter-dominated phase to the late-times accelerated one and in [27] the authors consider 
cosmological scenarios based on f(R,T) theory and the function f(R,T) is numerically reconstructed 
from holographic dark energy. 

In this paper, we consider the special case in which the function f(R, T) is the usual Einstein- Hilbert 
term plus a correction g(T), i.e. f(R, T) = R+ 2g(T). A suitable expression is assumed for the Hubble 
parameter, which may provide the four FTFS. For some values of the parameter a, differential equations 
of g(T) are established and solved. Another aspect which we study are quantum effects due to conformal 
anomaly near the singularities. We observe that conformal anomaly cannot remove the sudden singularity 
or the type IV singularity. However, for some values of the parameter a, the big rip and the type III 
singularity may be avoided. Another important result that we obtain is that, even without taking into 
account quantum effects, the big rip and the type III singularity may be avoided thanks to the contribution 
T of the modified theory. 

The paper is organized as follows. In Sec. 2, we present the general formulation of the theory 
addressing the special case f(R,T) = R + 2g(T) and obtaining the models that reproduce each type of 
FTFS. In Sec. 3, we investigate quantum effects due to conformal anomaly We present our conclusions 
and perspectives in Sec. 4. 



2 f(R, T)-gravity 

2.1 The model 

We assume a modification of general relativity in which the Ricci scalar R is replaced by an arbitrary 
function f(R,T). The total action reads 

S=^Jd i x^lf(R 7 T) + £}, (1) 

where C is the matter Lagrangian density and T = g^T^ is the trace of the matter energy-momentum 
tensor T^ v which is defined as 

_ 2 S(V-gC) 

Varying the action (TTJ) with respect to the metric, one obtains the gravitational field equations 



f R (R, T)R^ - -f(R, T) 9llv + ( 9llv n - V M V„) f R (R, T) = T^ - f T (R, T)T^ - f T (R, T)0^ , (3) 



where f R and /y denote the derivatives of / with respect to R and T, respectively; V^ is the covariant 
derivative and M! , is defined by 

o - a a0 5Ta P - 2T +o £ 2a afi - Lm U) 

^ ~ Sg^ v ~ £1 nv-^9iivl-m *9 dg^Qcc/3 ■ W 

We adopt the flat Friedmann-Lemaitre- Robertson- Walker (FLRW) metric ds 2 = dt 2 — a 2 (i)dx 2 , as the 
space-time geometry, and we assume the matter content of the universe to be a perfect fluid, whose the 
stress tensor is given by 

Tp, = (p + p)u^u v - pg^u , (5) 

where u^ is the four- velocity and satisfies the condition u^u^ = 1; p and p are, respectively, the energy 
density and the pressure of ordinary matter. The matter Lagrangian density can be taken as C m = —p, 
which implies 0^„ = —2T >1V — pg^v Thus, Eq. ([3]) becomes 

f R (R, T)R„ V - lf(R, T) 9flv + (g^D - V M V„) f R (R, T) = T M „ + f T (R, T)T^ + p 9flv f T (R, T) . (6) 

Our ansatz for the function / is the following: f(R, T) = R + 2g(T), where g(T) is an arbitrary function 
of T. Then, Eq. (JB) becomes 

R^ - \R9h» = T„u + 2g T (T)T^ + [2pg T {T) + g(T)] g^ . (7) 

The components 00 and ii of Eq. ([7]) read 

m 2 = p eff = p + 2(p + p)g T (T)+g(T), (8) 

-2H - 3H 2 = p eff = p - g{T) . (9) 



Here, p e // and p e // are the effective energy density and pressure, respectiveiy. Combining Eqs. (jHJ and 
dU), we get the following equations 



2H = p(l + uj)[l + 2g T {T)} , (10) 



, « 



6- = 2 [p + S (T)] - p (1 + w) [2ffr(T) + 3] , (11) 

a 

where we have defined u = p/p- 

Note that, in order to provide an accelerated expanding universe, the function g{T) has to satisfy the 
condition 

2g(T) - 2p (1 + u) g T (T) > p (1 + 3w) . (12) 

2.2 Finite-time singularities in R + 2g(T) gravity 

We are interested in the modified f(R, T) gravity models that produce some types of finite-time future 
singularities (see [551 US]): given by the Hubble parameter 

H = h(t s - t)- a , (13) 

where h and t s are positive constants and t < t s . This expression is chosen for guaranteeing an expanding 
universe. Due to fact that H or some of its derivatives (and therefore the curvature) could become singular 
when t is close to t s , the parameter a could be either positive or negative. If a > 0, then H diverges as the 
singularity time is approached. The same expression of the Hubble parameter provides the acceleration 
of the universe, since the strong energy condition is violated. On the other hand, when a < 0, a relevant 
positive constant H s may be introduced and interpreted as the Hubble parameter at the singularity time, 
such that H = h(t s — t)~ a + H s . In this case, some derivatives of H (and therefore the curvature) become 
singular. Note the particular case a = 0, which corresponds to the de Sitter space (which we do not treat 
in this paper). 

Depending on the value of a, different expressions can be found for the scale factor. 

If a = 1, the scale factor reads 



a(t) = a{t s - t)~ h , (14) 



where a is an integration constant. 
If a > and a ^ 1, we get 



a(t) = aexp 



h(t a - 1) 



l-ct 



a — 1 
Finally, if a < and H s > 0, the scale factor behaves as 

h(t s -ty 



a(t) = aexp < — (t s — t) 



H H 



1 



(15) 



(16) 



• 



• 



• 



The finite future time singularities can be classified as follows 

• Type I (big rip): for t — > t s , a — > oo, p e ff and \p e ff\ — > oo at t — t s . This corresponds to a > 1. 
Type II (sudden): for t — > t s , a — > a s , p e ff -^ Ps and |p e //| —* oo. It corresponds to —1 < a < 0. 
Type III: for t — > t s , a — > a s , p e ff —> oo and |p e //| - ► oo. This corresponds to < a < 1. 
Type IV: for t — > t s , a — > a s , p e ff — > 0, p e ff — > and higher derivatives of H diverge. This 
corresponds to the case a < — 1, but with a not an integer number. 

We investigate f(R,T) = R + 2g(T) gravity models which generate the above types of finite-time 
singularities. This corresponds to determine which form of g(T) gives rise to each type of singularity. 

We introduce an arbitrary auxiliary field which should not be confused with the scalar field on 
which the matter action depends. Write g(T) as 

g(T) = P{<P)T + Q{<t>) . (17) 

By using Eq. (|T7|) in Eq. (JXJ) with the ansatz f(R,T) = R + 2g(R) and varying the action with respect 
to the field <fi, one obtains 

P'(£)T + Q'(<£) = , (18) 

where the prime denotes the derivative with respect to the field <f>. Here, it is considered that the action 
is unchanged under the variation of the auxiliary field. Then, the derivative of the function g(T) can be 
written as 

g T (T) = P(c/>) + [P'(ct>)T + Q' (</>)] ^ . (19) 

Making use of Eq. (fT5)l , one gets 

g T (T) = P(0) . (20) 

From Eq. (jlOl) . one obtains 

P&) = ~\ ~ -rr—, , (21) 

2 p{l+Lj) 
and from Eqs. ©, dTTJ) and ([2"DJ). one finds 

Q(4>)=3H 2 -p-p(3-u)P(cf>). (22) 

Note that, since h is always positive for guaranteeing an expanding universe, an analysis can be done on 
the basis of ui. Then, one can observe two cases: u> > —1 and uj < — 1. The latter (a; < —1) is particularly 
interesting because even if g(T) vanishes, singularities may appear. 

I. The first case: u < — 1 



2.2.1 Treating the case a(t) = a(t s - i)~ /l 

With this scale factor, the energy density can be written as 



Substituting Eq. (|22|) into Eq. (J2T1) . one obtains 

1 ha^ 1+ ^ 

Now, injecting Eqs. (gSJ) and dH} into Eq. $n§, we get 
2 , ft(3-w) 



(t. - t)- 2 - 3hil+ ^ 



Q{4>) 



3h z 



(t s ty 2 + \ (1 - U ) P0 a~ 3 ^ (t, - t) 3/l(1+w) 



(23) 



(24) 



(25) 



1 + CJ 

Since one can redefine <fr properly, we may choose it to be the time coordinate, <fi = t. 

For u < — 1, being h > 0, one has that 3/i(l + w) is always negative. In such a situation, three 
conditions can be considered: 

• -2 < 3/i(l+w) < , 

• 3h(l + w) = -2 , 

• 3/i(l + w) < -2 . 

Let us now analyse each expression of the scale factor, finding the corresponding g(T). 

A The first condition: 2 < 3/i(l + w) < 

In this case, as the singularity time is approached (t — >• t s ), one has 



and 



Pit) ~ 



Q(t) 



(t 3 - t) 



p Q (l+ui) 

2 /l(3 — Ul 



-2-3/1(1+01) 



3/i 2 



(t. - t)" 



(26) 



(27) 



1 + w 
On the other hand, we can write the trace of the energy-momentum tensor in terms of the energy density 

as 



T = p-Zp 
= p(l-3u) 

from which one has p = T/(l — 3w), and then 

t s -t = ai [p (l - 3w)] _3ft < ll +") T 3 mi+-) . 

Now, substituting Eq. (f2"9l into Eqs. (|2"o) and (I2T1) . one obtains 



Q(i) -► g(T) 



a h 



1 +LO 

2 , ft(3-w) 



3/1 



1+UJ 



[p (l -3a;)] 3h d+-) T 3h d+-) 



(28) 

(29) 

(30) 
(31) 



Using Eqs. ^TTJ), ((30j) and ((3TJ), we get 

g(T) = ha-^ (3h + 2)[p a (l-3uj)}^^T^^y . (32) 

Thus, when the condition —2 < 3/i(l + w) < is satisfied, the big rip can appear with the model 

f{R,T) = R + 2hcT% (3ft + 2) [p (l - 3u)}^+^ TWtt=j . (33) 

B The second condition: 3/i(l + uj) = —2 

In this case, as t — > t s (considering t s — t + e, with e << 1), using Eq. (|29|) . Eqs. ([24)) and (|25|) become 

1 3/i 2 a~£ 



P(t)^P(T) = -- 



2po 



Q(i) -> Q(T) 



l + 3h 



1 h 2 a~^ 
3 po 



r 



(34) 

(35) 



2(l + 2/i) 

Now, substituting Eqs. (l34|) and (J35I) into Eq. (|T7|) . (/(T) is found and the corresponding f(R,T) model 
reads 



f(R,T) = R 



2 + 3/i 



3ti 2 a~t 



1 



Po 



T 



(36) 



3(1 + 2h) 

Hence, when the condition 3h(l + us) = — 2 is satisfied, the big rip can occur with the model of Eq. ([36 
C The third condition: 3/i(l + w) < -2 

In this case, as the singularity time is approached, using Eq. (f2"!)|) . Eqs. (IMt and (j2"5)l become 

1 



P(£) -> P(T) = 



2 ' 



Q(t)^ = i(l- W ) Po a- 3 ( 1+ -)(/ s -i) 3 ' ,(1+w) 



Q(T) 



1-w 



■T . 



(37) 



(38) 



2(1 - 3w) 

Hence, substituting Eqs. (|37p and (I3"51) into Eq. (fTTj) . g(T) can be obtained and the corresponding f(R, T) 
model reads 

2w 



f{R,T) = R- 



■T . 



1-3cj 
This model also allows the occurrence of the big rip with h satisfying the condition 3/i(l + cj) < —2 



(39) 



h(t s -ty 



2.2.2 Treating the case a(t) = a exp 

In this case, the appearance of both the big rip (a > 1) and type III singularity (0 < a < 1) is 
possible. The expressions ([21]) and (|22|). in this case, read 



P(t) 



1 aha 3{l+ ^ (t s - t) 



-ct-l 



/9q(1 + w) 



■ exp 



3/i(l+ u)(t s -t) 1 



a- 1 



Q(t) = 3h 2 (t s ~ty 2a + 



ah(3 — Ld)(t s — t)~ 



+ -{l-u})p a 



-3(l+w) 



exp 



1 + w 

-3/i(l + w)(i s -t) 1 - Q 
a-1 



(40) 



(41) 



Making use of Eq. (|28|) . one gets 

i, -t = ln~ 



[p (l - 3w)] 3 M 1 +")a i ^ T 3 ^ 1 ""' 



(42) 



• Type III singularity, < a < 1 

As the singularity time is approached, we get 

-aha^+^its-t)-*- 1 



Pit) 



Q(t) = 



/9o(l + w) 

*/i(3 - w) (t s - t) 
1 + ui 



-a-l 



(43) 
(44) 



Substituting Eq. g2|) into Eqs. (gSJ and (|J4)l. and using Eq. ([TT|). one can obtain g(T), and the f(R,T) 
model that leads to the type III singularity reads 



f(R,T) = R 



2a/ia 3(1+w) 
p (l +w) 



go(3--a0 

«3(l+w) 



-T 



ln^ 



a-l 1-c, 1-a 

[po(l - 3w)] 3 mi+-> a~^~ T 3fc ( 1 +") 



• Big rip, a > 1 

In this case, as the singularity time is approached, Eqs. (J40I) and (l4Tj) become 

afco 3 ( 1+fc 'H*«-*)" a ~ 1 



P(t) 



Po(l+w) 



1 



Q(*) ~ ^( 1 -w) y o a 



-3(l+w) 



exp 



-3/i(l + w)(t s -t) ] 



On the other hand, from Eq. (|4"2"j) we can write 

-3/i(l + w)(i s -t) 1 -" 



exp 



:- 1 



a — 1 



=3(l+o;) j. 



from which, we get 

P(t) -► P(T) = 

^-O^ - 2(1 _3a,) 
and <?(T) can be obtained through Eq. (|T7| . leading to the model 



ate 3(i+ w ) a+1 
1-w 



po(l-3w) 

Q-l 1-q 1-Q 

[po(l - 3a;)] 3fe ci+") a - ? 1- T 3fc ( 1 +") 



■r , 



l-w 2a/ia 3 ( 1+ ") . «+!_! 



In" 



I 1 — 3w po{l + u>) 
which can allow the occurrence of the big rip. 

2.2.3 Treating the case a(t) — a exp \—(t s — t) 

Here, the expressions (|2"TT) and (l2"2"j) take new forms 
p(t) _ 1 aha 3( - 1+ "\t s - 1)-"- 1 

Q(t) 



[p (l - 3w)] 3 Mi+") a 1 ^ T 3fe < 1 +"' 



T 



#o 



h(t s -t)' 



a-l 



3[h{t s -t)- a +H s ] 



exp i -3(1 + u)(t.-t) 
po(l + w) t 

2 a/i(3-w)(t s -i)"" -1 



#« 



h(t s -t)- a 



a — 1 



1 +CJ 



+ 2/ 9 o^ 3(1+ " )ex P1 3 ( 1+w )(^- f ) 



ff, 



Kt s -ty 

a-1 



(45) 



(46) 
(47) 

(48) 

(49) 
(50) 

(51) 



(52) 



(53) 



• Sudden singularity, — 1 < a < 0. 

As the singularity time is approached, Eqs. ([52]) and (f53|) can be written as 

P(t) 

Q(t) 



Po(l+«) 
ah(3-w)(t e -t)- l - a 



1+CJ 

On the other hand, the trace can be put in the following form: 



T = p (l - 3w)a -3(1+ ' j) exp \ 3(1 + oj)(t s - t) 
which, near the singularity, leads to 



H s - 



h(t s - ty 

a — 1 



ts — t = In 



\p0(l-3w)] 3H,(l+»)aH,J 1 3H s (l+.) 



(54) 
(55) 

(56) 

(57) 



Substituting Eq. (1571) into Eqs. (|54"|) and (|55|) . and using Eq. (|T7)) . g(T) is obtained and the model 
corresponding to the sudden singularity reads 



f(R,T) = R + 4ah\n 



[po(l-3w)j 3H 3 (i+")a 7 ^r 3 H s ( 1 +-) 



(58) 



• Type IV singularity, a < — 1. 

In this case, close to the singularity, Eqs. (|52|) and (|53| become 



P(t) 



1 
2 ' 



Q(<) ~ 3if s 2 + -(l-u;)poa- 3(1+w) 



(59) 
(60) 



Here, as the singularity time is approached, the trace of the energy-momentum tensor takes on the simple 
expression 



T = p a- 3( - 1+ ^(l-3Lu) 



(61) 



from which, using Eq. ([T7|). g(T) is easily obtained and then, the f(R,T) models which allows the 
occurrence of the type IV singularity reads 



f(R,T) = R 



2w 



l-3w 



T + QHi 



(62) 



II. Commenting the case u > — 1 

• Considering the scale factor of Eq. ()14|) 

In this case, the quantity 3h(l+ui) is always positive, such that — 2 — 3h(l + u) < 0. This leads to the 
same expressions as that in Eqs. (|2"6"|) and (j2"7) for P(t) and Q(t), respectively. Then, the corresponding 



f(R,T) model is already given in Eq. (|33|) . Therefore, the big rip can also appear, even if u> > —1. It 
is important to note that if g(T) vanishes, the fluid corresponds to quintessence and then, the type III 
cannot occur. Our result is clearly the effect of the presence of a part depending on the trace of the 
energy-momentum tensor in the gravitational action. 

• Considering the scale factor of Eq. (|15p 

In this case, for < a < 1, corresponding to the appearance of the type III singularity, the expressions 
for P(t) and Q(t) are the same as in Eqs. (|43)) and (|44|) and then the corresponding f(R, T) model is the 
one in Eq. ((45]). 

For the case a > 1, corresponding to the appearance of the big rip, the situation is quite different. As 
the singularity time is approached, the expressions (|40|) and (|4Tj) become Eqs. ([54]) and (|55|) . respectively, 
and the corresponding f(R, T) model is the same as in Eq. ([58]) . 



• Considering the scale factor of Eq. (|16|) 

For — 1 < a < 0, the case in which the sudden singularity may appear, Eqs. (|52[) and (|53|) take the 
form (|54j) and (|55j) and consequently, the corresponding model is given in Eq. f|58[) . 

For a < — 1, the case in which the type IV singularity may appear, the expressions (|52[) and ()53[) 
become (|59[) and (|60|) , respectively, and the corresponding model is given in Eq. (|62|) . 



3 Quantum effects near finite time singularity 

We may check the avoidance of finite time singularities by taking into account quantum effects. In 
this work we address conformal anomaly. Near the future finite-time singularity (t — > t s ) the curvature 
diverges. Due to their dependence on the curvature, quantum effects coming from conformal anomaly 
also become important. In such a situation, all classical considerations have to be revised and any claim 
about the appearance of future finite-time singularity cannot be justified without an account of quantum 
effects. One may incorporate the massless quantum effects by taking into account the conformal anomaly 
contribution as a backreaction near the singularity. The conformal anomaly Ta has the following well- 
known expression |32j 

T A = b f F + ^DRj + b'G + b"UR , (63) 

where F is the square of the four-dimensional Weyl tensor and G is Gauss-Bonnet invariant, given by 

F = ^R 2 -2R llv R» v + R liVpa R^r° , (64) 

G = R 2 - AR^R^ + R^Rf""" . (65) 



10 



Explicitly, if there are N scalars, N\/2 spinors, N\ vector fields, N 2 gravitons and Nhd higher derivative 
conformal scalars, b and b' are expressed as 



N + 6N 1/2 + 12iVi + 6IUV2 - 8N HD 
1920tt 2 



6 = TT^Za . ( 66 ) 



N + HJV 1/2 + 627Vi + 1411JV 2 - 287V ffD 
5760tt 2 



& ' = — ^ . ; m , • - - ■" l( r,i 

whereas 6" is an arbitrary constant whose value can be shifted by the finite renormalization of the local 
counter-term R 2 . 

Now, by taking into account quantum contribution from the conformal anomaly, Eq. ([7]) is modified 
and its trace can be written as 



-R=[l + 2g T (T)}+8pg T {T) + 4g{T) + T A . 
From now on, we assume that the ordinary fluid is pressureless, i.e. p — 0. Then, Eq. ([55]) becomes 

-R=[l + 2g p (p)]+4g(p)+T A . (69) 

For the FLRW universe, it is clear that 

F = , G = 24 (ilH 2 + H 4 ) . (70) 



Let us focus our attention on the case 2b + 36" = and set V = 1. Then, one gets Ta — G. We may 
analyse quantum effects coming from conformal anomaly near each type of singularity. 

3.1 Conformal anomaly near the big rip 

In this case, one can observe two conditions on the parameter a, that is, a = 1 and a > 1. 

• The case a = 1 

Here, the curvature behaves as R oc (t s — t)~ 2 , the trace as T ex (t s — ^) 3 ' l ( 1 +") j arLC l T A ^ (t s — t)~ A . 

Here, we may distinguish the three conditions previously mentioned; —2 < 3/i(l +uj) < 0, 3h(l +u>) = 
-2 and 3/i(l+ w) < -2. 

For -2 < 3/i(l+a;) < 0, the corresponding expression of g(T) is Eq. (J32J) , and then, g(T) oc Tg T {T) oc 
(t s — i)~ 2 . It appears that, as the singularity time is approached, the conformal anomaly term dominates 
other others terms. Hence, the big rip may be avoided. 

For 3/i(l + uj) — —2, the corresponding expression for g(T) is that coming from Eq. (1361) . and then 
g(T) oc Tgx{T) oc (t s — t)~ 2 . Also, in this case, the conformal anomaly term dominates and the big rip 
may also be avoided. 

For 3/i(l + to) < —2, the corresponding expression of g{T) is that coming from Eq. (|55|) . and g(T) oc 
TgxiT) oc (t s — i) 3 '»( 1 + w ). Here, it can be observed that, even if the conformal anomaly is not taken into 
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account the big rip may be avoided from the contribution of the matter in the gravitational part of the 
action. This can be viewed as a reason for taking into account the trace T term in the gravitational part 
since, in this case, the big rip can be avoided without the need of quantum effect. 

• The case a > 1 

Here, as the singularity time is approached, the curvature behaves as R oc (t s — t)~ 2a , the conformal 
anomaly as (t s — t)~ 4a , and the ordinary trace asToc e^ ts ~^ ^ " . The corresponding expression of 
g(T) is that coming from Eq. (f5Tj) and one gets, g(T) oc Tgx{T) oc (t s — t)~^ L e^ ts ~ t ^ " , that one 
can assimilate to (t s — t)~ 2N , where N — N(u,u) is a positive function of a and ui. This assumption is 
possible since, for a > 1, — - < 0, and also, as we are dealing with the case u < — 1, 3/i(l + ui) < 0. 
It is clear here that we also have —4a < —2a, and then, the conformal anomaly term dominates the 
curvature one. One can conclude that for a > 1, quantum effects from conformal anomaly may allow the 
avoidance of the big rip. Note also that, if for some values of a > 1 and ui < — 1, N > 2, then, we regain 
the situation on which f(R,T) itself lead to the avoidance of the big rip, from the ordinary trace term 
in the gravitational part of the action, without the need of quantum effects. 

3.2 Conformal anomaly near the sudden singularity 

In this case, —1 < a < 0, and as the singularity is approached, the curvature behaves as R oc 
(t s — t) _Q_1 , the ordinary trace as T oc e^ 3- *^ " and the conformal anomaly as Ta oc (t s — t) _3QI_1 . 

The corresponding expression of g(T) is that coming from Eq. (1581) . Here, one gets g{T) oc Tgx(T) oc 
(t s — t)~ a ~ l . In this case, it appears that the avoidance of the sudden singularity must come from the 
quantum effect. In this way, one can observe that for —1 < a < —1/3, — 3a— 1 > while — a— 1 < 0. Then, 
the sudden singularity is robust again quantum effect coming from conformal anomaly. For —1/3 < a < 0, 
— a — 1 < — 3a— 1<0. Also in this situation, the sudden singularity cannot be avoided from conformal 
anomaly contribution. 

3.3 Conformal anomaly near the singularity of type III 

This type of singularity occurs for < a < 1, and around the singularity, the curvature behaves as 
R oc (t s — i) - " -1 , the ordinary trace as T oc e(* s- *) " and the conformal anomaly as Ta oc (t s — 1) — 3a— 1 . 
The corresponding expression of g(T) is that coming from Eq. (|45|) . and we gets g(T) oc Tgx{T) oc 
e'* s_t - ) . As < a < 1, one gets —3a — 1 < —a — 1 < 0. Then, the singularity of type III may be 

avoided from conformal anomaly effects. 
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3.4 Conformal anomaly near the singularity of type IV 

The condition of occurrence of this type of singularity is a < — 1, and this case the curvature is 
not singular, but form it first to higher derivatives diverge as the singularity time is approached. The 
conformal anomaly term behaves as (t s — i) _3Q_1 . Note that in this case, we always get —3a — 1 > 0. In 
such a situation the singularity is always robust against the conformal anomaly and cannot be avoided. 

Another important point to be observed is the case where u > — 1. The special aspect in this case is 
that the quantity 3h(l +u>) is always positive, and then, any singular behaviour cannot be observed from 
the ordinary trace contribution. Hence, the unique possibility for avoiding the singularities must come 
from quantum effects. 

4 Conclusion 

We considered the modified f(R, T) theory of gravity where R is the curvature scalar and T the trace of 
the energy-momentum tensor. We focused our attention in reconstructing the models of this theory, that 
may lead to finite-time future singularities. We assumed a suitable expression for the Hubble parameter 
in order to provide the expansion of the universe and realise the four finite-time future singularities, by 
adjusting an input parameter a. A special expression f(R, T) = R+2g(T) is assumed, viewed as a T-term 
correction to the Einstein-Hilbcrt term, R. From this, differential equation for g(T) are obtained through 
the equations of motion, for some values of the parameter a. In each case, the differential equation is 
solved leading to the f{R,T) model that leads to finite-time future singularity. 

As the curvature diverges at singularity time, quantum effect may be introduced. In this way, con- 
formal anomaly is introduced and its effects at singularity time in analysed. We observed that conformal 
anomaly cannot remove the sudden singularity nor the type IV singularity. However, for some values of 
the parameter a, the big rip and the type III singularity may be avoided. Another important result that 
we find is that, even without the contribution of quantum effect, the big rip and the type III singularity 
may be avoided from the contribution of the ordinary trace terms, for some values of the parameter a. 
This is a reason for considering f(R, T) in the study of some scenarios of evolution of the universe. 

In the case of type II (sudden singularities) and type IV singularities, we observed that neither 
conformal anomaly or the terms coming from the classical trace allow the avoidance of these singularities. 
Still, in order to check the possible avoidance of type II and IV singularities, the viscosity of the fluid 
characterising the matter content may be introduced. We leave this as a future work. 

Acknowledgement: M. J. S. Houndjo thanks IMSP (Benin) and UAC (Benin) for the hospitality 
during the elaboration of this work. 
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